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An abundance of invariant polynomials satisfying the 

Riemann hypothesis 

Koji Chinen* 
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Abstract 

In 1999, Iwan Duursma defined the zeta function for a linear code as a generating 
function of its Hamming weight enumerator. It can also be defined for other homoge- 
neous polynomials not corresponding to existing codes. If the homogeneous polynomial 
is invariant under the Mac Williams transform, then its zeta function satisfies a functional 
equation and we can formulate an analogue of the Riemann hypothesis. As far as exist- 
ing codes are concerned, the Riemann hypothesis is believed to be closely related to the 
extremal property. 

In this article, we show there are abundant polynomials invariant by the Mac Williams 
transform which satisfy the Riemann hypothesis. The proof is carried out by explicit 
construction of such polynomials. To prove the Riemann hypothesis for a certain class of 
invariant polynomials, we establish an analogue of the Enestrom-Kakeya theorem. 

Key Words: Zeta function for codes; Riemann hypothesis; Perfect code; Enestrom-Kakeya 
theorem; reciprocal equation; Invariant polynomial ring. 
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: 1 Introduction 

Let p be a prime, q = p r for some positive integer r and we denote the finite field with q 
elements by F q . Let C be an [n, k, <i]-code over F q with the Hamming weight enumerator 
Wc(x,y). Duursma [I] defined the zeta function for C as a generating function of Wc(x, y). 
Then the author [2] considered the case of so-called "formal weight enumerators", noticing 
that Duursma's definition can be extended for other homogeneous polynomials than the weight 
enumerators of actual codes. Taking these into account, we start from the following definition: 



Definition 1.1 For any q G N (q > 2) and any homogeneous polynomial of the form 

n 

W(x,y) = x n + J2 Ax n -y (A4 G C, A d ^ 0) (1.1) 

i=d 
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there exists a unique polynomial P{T) G C[T] of degree at most n — d such that 



P(T) 



( y (l-T) + xT) n = --.+ 



W{x,y)-x n T 
q-1 



m—d 



+ ••■• 



(1.2) 



(1-T)(1 



qT) 



We call P{T) and Z(T) = P(T)/(1 — T)(l — qT) the zeta polynomial and the zeta function of 
W(x, y), respectively. 

For the proof of existence and uniqueness of P(T), see Appendix. If W(x,y) = Wc(x,y) for 
some linear code C, then we take q in the above definition as JjF 9 , but if W(x, y) is not related 
to an existing code, then q must be chosen suitably according to what meaning W(x, y) has. 

In the case W(x, y) = Wc(x,y), the zeta polynomial P(T) for Wc(x,y) is of particular 
interest when C is self-dual, because it has the functional equation 



(g = n/2 + 1 — d, see [51 p.59]), which is a result of the fact that Wc(x,y) is invariant by the 
MacWilliams transform 



where we define f a (x, y) = f(ax + by, cx + dy) for f(x, y) G C[x, y] and a linear transformation 



The functional equation ( 11. 3D is the same as that of zeta functions of algebraic curves, so 
we can formulate the Riemann hypothesis (see Duursma [HI Definition 4.1]). Even if W(x,y) 
does not correspond to an actual code, we can formulate the Riemann hypothesis in the same 
way provided that W trq (x,y) = W(x,y) because it is this invariance that yields (11.31) : 

Definition 1.2 The code C (or the invariant polynomial W(x,y)) satisfies the Riemann hy- 
pothesis if all the zeros of P(T) have the same absolute value 1/y/q. 

Duursma deduces various interesting properties of P(T) and discusses their possible applications 
to the coding theory (see [5J EJ [7]). 

Finding an equivalent condition for the Riemann hypothesis above seems still an open 
problem, but Duursma asks the following ([6, Open Problem 4.2]): 

Problem 1.3 Prove or disprove that all extremal weight enumerators satisfy the Riemann 
hypothesis. 

A self-dual code C is called extremal if it has the largest possible minimum distance (see Pless 
[TT| p. 139]). There are 4 well-known sequences of extremal self-dual codes (Types I, II, III and 
IV, see Conway-Sloane [3]). The extremal code is also characterized by its weight enumerator 
Wc(x, y): the code C is extremal if d of Wc{x, y) in ( 11.11) is the largest among all the self-dual 
weight enumerators of degree n over F q . Using this, the extremal property is straightfowardly 
extended to the case of some more general invariant polynomials. Duursma proved that all 
extremal Type IV codes satisfied the Riemann hypothesis ([7]). Thus, as far as the existing 




(1.3) 




(1.4) 
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codes are concerned, we may expect that the Riemann hypothesis reflects one of the abilities 
of the code, the extremal property. 

In [2] , the author extended the consideration to the case of the formal weight enumerators. 
A formal weight enumerator W(x, y) resembles the weight enumerator of a Type II code, but is 
distinguished from it by the property W a2 (x,y) = —W(x,y) (see [21 Definition 1.4]). The zeta 
polynomial P{T) of W{x,y) satisfies P(T) = ~P{l/2T)2 9 T 29 (g = n/2 + 1 - d) and we can 
formulate the Riemann hypothesis in the same way as in Definition 11.21 setting q = 2. In [21 
Section 3], we observed that the extremal property might yield the Riemann hypothesis also in 
the case of the formal weight enumerators. 

The purpose of the present article is to extend the consideration to all the polynomials 
which are invariant by the MacWilliams transform a q . Such polynomials form an invariant 
polynomial ring 

C[x,y} Gq = C[x+yq-l)y,y(x-y)\ (1.5) 

where G q = (a q ) (see MacWilliams- Sloane [91 p. 605, Theorem 5]). As a problem of invariant 
polynomials, we can remove the structure of linear codes and allow q to be any positive integer 
such that q > 2. We try to find as many polynomials as possible in C[x, y] Gq which satisfy the 
Riemann hypothesis. The results imply that the Riemann hypothesis is not always relevant to 
the extremal property in the ring C[x,y] Gq . The first result is the following: 

Theorem 1.4 For any q > 2 and any n, d such that 2 < d < there exists a o q -invariant 
polynomial of the form U.l\) which satisfies the Riemann hypothesis. 

Note that the restriction 2 < d comes from the original Duursma theory. We also note that 
the number d in (11.11) must satisfy d < | + 1 in C[x, y] Gq . In cases where equality holds, 
the polynomial becomes MDS and the zeta polynomial is a constant (see Section 3). Thus 
by the condition 2 < d < almost all possible pairs of n and d are covered and it shows 
that polynomials satisfying the Riemann hypothesis are widely and abundantly distributed in 
C[x,y} Gq . 

The notion of extremal polynomial is defined only in terms of n and d, but Theorem 11.41 
implies that, at least in the ring C[x, y] Gq , the condition for the Riemann hypothesis is not 
determined by n and d only. Thus Theorem 11.41 shows us another aspect of the zeta functions 
for invariant polynomials. 

Theorem II .41 is proved by explicit construction of the invariant polynomials with the desired 
property. This is done by using the weight enumerators of codes which are not self-dual. If 
C is not self-dual, its weight enumarator Wc(x,y) does not satisfy W^^x^) = Wc(x,y), but 
combining Wc(x,y) and W c ±- (x, y) , we can easily get an invariant expression Wc(x,y) and its 
zeta polynomial Pc{T) (see Section 2). Theorem 11.41 is the result of the case where C is an 
MDS code (see Section 3). 

Such a way of constructing invariant polynomials can be applied to any linear code which is 
not self-dual and leads to further exploration. The rest of the paper is devoted to the analysis 
of two other special classes of codes, the general Hamming codes and the Golay codes (not 
self-dual). These codes, along with certain MDS codes form an important class of good codes, 
the perfect codes (see Pless [HI p.21]): 

Definition 1.5 A code C C F q n of minimum distance d is called perfect if all the vectors in 
F g n are contained in a ball of radius [(d— 1)/2] about the codewords, where [x] means the largest 
integer not greater than x. 
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The nontrivial linear perfect codes are completely determined ( [TTI Section 2.2] or [91 Section 
6.10]): 

(i) The general Hamming [(q r — l)/(g — 1) = n, n — r, 3] codes over F g , 

(ii) The binary [23, 12, 7] and the ternary [11, 6, 5] Golay codes. 

We also have trivial perfect codes: the whole space and a binary repetition code of odd length. 
The latter has the parameter [n, l,n], being MDS and dealt with in Theorem 11.41 As to the 
general Hamming codes, they become MDS when r = 2, so it follows that this case is also 
treated in Theorem 11.41 

We can find again infinitely many polynomials in C[x, y] Gq satisfying the Riemann hypoth- 
esis by constructing Wc(x,y) from the above class of codes: 

Theorem 1.6 Let C = Ham(r, q) be the Hamming [(q r — l)/(q — l) = n, n — r, 3] code over F q . 
If r > 3 and q > 4, then the invariant polynomial Wc(x, y) in C[x, y] Gq satisfies the Riemann 
hypothesis. 

To prove Theorem 11.61 we deduce a certain function theoretical result concerning the distribu- 
tion of the zeros of a self-reciprocal polynomial: 

Theorem 1.7 If f(T) = a + ai T + ■ ■ ■ + a k T k + a k T m ~ k + a k ^T m ~ k+1 + ■■■ + a T m (m>2k) 
satisfies > a\ > ■ ■ ■ > a k > 0, then all the roots of f(T) lie on the unit circle. 

This is, so to speak, a self-reciprocal analogue of the famous Enestrom-Kakeya theorem (see 
Theorem 1 5. II) . Because of the technical difficulties, Theorem 1 1 . 61 remains unproved when q = 2, 3 
and r > 3, but numerical experiments imply that the Riemann hypothesis seems to be true in 
these cases. 

For the Golay codes, we have the following: 

Theorem 1.8 Let C be the binary [23,12,7] or the ternary [11,6,5] Golay code. Then the 
invariant polynomial Wc{x,y) satisfies the Riemann hypothesis. 

Thus except for the binary and ternary general Hamming codes, we can prove that the 
invariant polynomials Wa(x,y) from the perfect codes satisfy the Riemann hypothesis. 

The rest of the article is organized as follows. In Section 2, we construct an invariant 
polynomial Wc(x, y) from the weight enumerator of a code C (which is not always self-dual) 
and give an explicit form of its zeta polynomial Pc{T). In Section 3, we apply the results 
in Section 2 to the MDS code and prove Theorem 11.41 In Section 4, we determine the zeta 
polynomial Pc{T) when C = Ham(r, q), the general Hamming code when r > 3 and q > 2. 
Section 5 is devoted to an analogue of the Enestrom-Kakeya theorem. Here we use several results 
of the classical function theory. Using it, we prove the Riemann hypothesis for H'Hamfr.g) ( x i u) 
{t > 3, q > 4) in Section 6. In Section 7, we consider the case of the Golay codes, and prove 
Theorem 11.81 by a different method to Theorem 11.61 

We have been interested in the extremal property of the weight enumerators when consid- 
ering the Riemann hypothesis in the context of existing self-dual codes or a little larger class 
of invariant polynomials which have some connections to the coding theory, that is, the formal 
weight enumerators. But the results in this article show that it is not always the extremal prop- 
erty that yields the Riemann hypothesis in the "largest" ring C[x, y\ Gq . We can observe rather 



4 



pathological phenomena there. We are now in a position to seek some new structures which 
are larger than existing codes (but smaller than C[x, y] Gq ), in which the Riemann hypothesis 
indicates some distinguished properties of invariant polynomials. 

Acknowledgment. The author would like to express his sincere gratitude to Professor Leo 
Murata for an abundance of valuable advice and discussion. 



2 Invariant polynomials and their zeta functions from 
arbitrary linear codes 

Let C be a linear [n, k, d] code over F g and Wc{x, y) be its Hamming weight enumerator. 
Suppose the dual code C x has the parameter [n, n — k, d 1 -) and we assume d, d 1 - > 2. Com- 
bining Wc{x,y) and the dual weight enumerator Wc±{x,y), we can easily obtain an invariant 
expression Wc{x, y): 

Proposition 2.1 Let 

Wc(x,y) := i - l qk _ n/2 {W c (x, y) + q k - n ' 2 W c . (x, y)}. (2.1) 

Then we have Wc q (x,y) = Wc(x,y), i.e. Wc{x,y) G C[x,y] Gq . 
Proof. The proof is evident from the MacWilliams identity 

W£(x,y) = q k - n ' 2 W c ^x,y) or W%_(x,y) = q n ^- k W c (x,y) 

(see 0, p. 146, Theorem 13]). | 

Now we deduce the explicit form of the zeta polynomial Pc(T) of Wc(x,y). Let Pc(T) 
and P C ±(T) be the zeta polynomials of Wc(x, y) and W c j-(x, y), respectively. Our goal in this 
section is to prove the following: 

Theorem 2.2 The zeta polynomial Pc{T) of Wc{x,y) is given by 

P C (T) = - - gfc . n/2 |p c (T) + q ^- d P c T^ 2 - M | . (2.2) 



It satisfies deg Pq = 2g and the functional equation 

A? 

where g := n/2 — 1— min(<i, d- 1 ). 



P C (T) = P c ( 4f ) Q S T 2 ~ 9 (2.3) 



Proof. By Definition 11.11 We have 



p c(T) «v , W c (x,y)-x n rrn _ d 



(y(l - T) + xT) n = ■■■+ ° v T n ~ d + ■■■ (2.4) 



;i-T)(l-gT)^ ' ' q-1 
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and 

_^L_ fe(1 _ T) + xTT = ... + ^M-^ + ■ • • . (2.5) 
We suppose < rf ± . Then (E2D multiplied by q k-n/2 T d -d 

becomes 

■^g^ W - T) + *r)» = . . . + ^"<WM-^ + . . . . (2 . 6) 

We add (J23D and (l2~^]) . then divide it by 1 + q k ~ n / 2 . It gives 

{Pc(T) + q k '^P c 4T)T d± - d }/(l + q k -^) , ^,g)-^ , 

(l-T)(l-gT) 12/1 JJ + ^J -•••+ ? _ 1 +•••• 

Thus we have 

Pc ^ = i + lk-n/2 { P d T ) + q k - n/2 Pc4T)T d± - d } (2.7) 

by the existence and uniqueness of the zeta polynomial. The polynomial P C ±(T) can be sub- 
stituted by 



where 



g = n + i-k-d, (2.8) 
g ± = k + l-d ± . (2.9) 

These formulas come from the original Duursma theory (see Duursma [5.; p.59]). Hence we 
have 

P ^ = I + qk-n/2 { P c< T ) + q n/2+1 - d Pc Qp) T" +2 - 2d | . (2.10) 
When d > d,- 1 , similarly we have 

P C (T) = - T + d q l n/2 \P C {T) + q ^~ d P c T^} . (2.11) 

These two formulas give (12.21) . The functional equation (12. 3p is obtained in a similar manner 
to that of Duursma [HI P- 119]. As to degP^, first we note that 

deg P c = deg P c ± = g + g 1 = n + 2 - d - d 1 - (2.12) 

(see [SJ p.59]). By (12. 7p . we have degP^ = n + 2 — 2d = 2g when d < d- 1 . The case d > d L is 
similar. | 

Remark. When C 1 - = C, we can easily verify that Pc{T) = Pc(T). Thus we have extended 
Duursma's theory in such a way that the zeta functions for codes which are not self-dual have 
the functional equation. 
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3 The MDS codes 



We consider the case where C is an MDS code in the construction of Wc{x, y) in Section 2 and 
prove Theorem ll.4[ An [n, k, d] code C is called an MDS (maximal distance separable) code if 
d = n — k + 1 is satisfied, i.e., the equality holds in the Singleton bound d < n — k + 1. If C 
is MDS, then so is C 1 - and it has the parameter [n, n — k, n + 2 — d]. The weight enumerator 
Wc{x, y) of an MDS code C is determined only by n, d and q. It can be explicitly given in 
terms of binomial coefficients: 

Theorem 3.1 Let Wc(x, y) = ^2i=o \x n ~ % y % be the weight enumerator of an [n, k,d = n—k+1] 
MDS code C . Then we have 

^=(I)E(-1)'(3(9" +1 - 3 '-D- «2f> 
Proof. MacWilliams-Sloane [91 p. 320, Theorem 6]. | 

We allow A{ to be negative and q be arbitrary integer greater than one. Even in the case 
Wc(x,y) does not represent the weight distribution of an actual code, we are interested in the 
polynomial itself and often call it an "MDS polynomial". From now on we assume d,d ± > 2. 
What is crucial for our discussion is the following: 

Theorem 3.2 Let C be MDS. Then we have P C (T) = 1. 

Proof. See Duusrma |5, Proposition 1]. In cases where q is not a prime power, it is straight- 
forward. | 

Now we determine the range of d and n. Both C and C L are MDS, so we can assume d < d L 
without loss of generality. Since d 1 - = n + 2 — d, d < d 1 - is equivalent to 

n 

rf<- + l. (3.1) 

If equality holds in (13.11) . then g = and P C (T) is a constant (Wc(x, y) is an MDS polynomial 
in the ring C[x, y] Gq in this case. It can happen when n is even). We exclude this case and 
have d < (n + l)/2. Duursma's theory requires d, d x > 2, therefore d and n can assume the 
values with 

n + 1 , . 

2<d<^— . (3.2) 

Now let C be an (actual or virtual) MDS code. Then we have Pc{T) = Pc±(T) = 1 by Theorem 
13.21 The zeta polynomial Pc{T) of the invariant polynomial Wc{x, y) is given by Theorem 12.21 

as 

Pc{T) = 1+g \_ n/2 (l + ^2+1-^+2-2^ (33) 

We can easily see that all the roots of (13.31) lie on the circle \T\ = 1/y/q- From Theorem 13.11 
and (13.21) . Wc(x, y) is of the form x n + AdX n ~ d y d + ■ • ■ and Ad ^ 0. This completes the proof 
of Theorem ll.4[ 
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4 The general Hamming codes 

For r > 2 and a prime power q, the general Hamming [(q r — l)/(g — 1) = n, n — r, 3] code 
Ham(r, q) over F g is the dual code of an [n, r, g r_1 ] simplex code over F q (see Pless et al. [12j 
p. 316]). Therefore we have 

n-l (q-l)n + l 

^Ham(r,g)^(a;,y) = l" + (?-l)ra ' J * (4.1) 

= x n + (g r -l)a; n - 9r " 1 i/ <?r " 1 - 

In this section we assume r > 3, allow g to be any integer with q > 2 and determine explic- 
itly the zeta polynomial P r ^ q (T) := P-Ra,ra(r, q )iT) of the invariant polynomial Wua,m(r,q)(x,y) G 
C[x,y] Gq constructed in the manner of Section 2. For our purpose, it is easier to handle with 
W / Ham(r, g ) ± ( ::c ) u) than Wnara{r,q) (x, y) , so we fix the notation as follows: 

C = Ham(r, q)- 1 , C 1 - = Ham(r, q), 
q r — 1 

n = (the length of C and C -1 ), 

q-1 

d = g 7 '" 1 (the minimum distance of Ham(r, q)' 1 ). 

First we deduce the zeta polynomial Pc{T) = P Ham ( r (? )±(T). We use the notion of the normal- 
ized weight enumerator (see Duursma [3 Definition 2]): 

Definition 4.1 For a weight enumerator A(x,y) of the form (CQJ), the normalized weight enu- 
merator a(t) is defined by 



-d 



Q 

The following theorem gives the relation between A(x, y) and its zeta polynomial P(T): 

Theorem 4.2 (Duursma) The weight enumerator A(x,y) , its zeta polynomial P(T) and the 
normalized weight enumerator a(t) are related by 



P(T) 



T) d+1 = a (jZTf^ ( mod T n ~ d+l ). 



{\-T){i-qry 

Proof. See [3 Theorem 2]. | 

For our code C = Ham(r, q)- 1 , the normalized weight enumerator a(t) is quite simple: 
Lemma 4.3 Let a r>q (t) be the normalized weight enumerator o/Ham^g) 1 . Then 



a r , q (t)=n I ( n _ 1 



i.e., a r ,q{t) is a constant. 



The proof is easy from (14. lft . Using this lemma and Theorem 14.21 we can deduce the explicit 
form of Pc(T): 
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Proposition 4.4 For r > 3 and q>2, the zeta polynomial Pc{T) = Pna,m<r,q) x (T) is given by 




where n/ ( r " i) ■ 

Proof. Lemma [4.31 and Theorem 14.21 gives 

P C {T) = a(—) (1 7 T)(1 "f ) (modr^) 

= ^^^(modT^ 1 ). ( 4 - 3 ) 

We have 

deg Pc = n + 2 — d — 3 = n — d — 1 < n — d + 1 

(see (I2.12p ). so Pc{T) coincides with the power series expansion of N r> g(l — qT) /(l — T) d up to 
the term of T n ~ d ~ l . By the expansion (1 — T)~ d = ("'d-T 1 )^' we have 



1 - qT 
(1-T) d 




This formula gives the desired result. | 

Remark. In the formula f fl~4l) . - q( j ^_ 2 ) = holds if and only if j = n - d. Thus the 

term of T n_d really vanishes in (14.41) . 

The main theorem in this section is the following: 
Theorem 4.5 For r > 3 and q > 2, the zeta polynomial P rj(? (T) := PHam(r,,j) (T) is given by 

Pr.iT) = ^^mT) - qF 2 {T)\ 

where 

— ?. — — 3 

n-d-2 , . x n-4 / • , \ 

*n - e T7> i+W2ri +E C- >*• 

Remark. If r = 2, both Ham(r, g) and Ham(r, q) 1 - are MDS codes and are treated in Section 
3. 

Proof. Since d = q r ~ x > 3 if r > 3 and q > 2, we have from Theorem 12.21 

PrAT) = i + T g tl /2 {^(T) + g"/ 2+1 - d P c f-L) T" +2 - M } . (4.5) 
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The remaining task is to describe each term in (I4.5p explicitly. We have from Proposition 14.41 



T 



d-3 



1 + q r - n / 2 



Pc{T) 



Nr, q 



1 + q r - n / 2 



n—d—l 



N, 



r.q 



1 + q r ~ n / 2 



rpd-3 _j_ 

n— 4 



i=d-2 



j + d - 1 

d- 1 



i + 2 
d- 1 



j + d-2 
d-l 



rpd+j- 



i + l 
d-l 



by putting d + j — 3 = z. Next we have from Proposition 14.41 again that 



T 



d-3 



1 + q r - n / 2 



■ q n ' 2+l - d Pc { 4 ) T 



qT 



-m+2-2d 



N r , q q n / 2+1 - d 
1 + q r ~ n / 2 



N, 



r,q 



1 _|_ q r-n/2 
n—d—l 

+ E 



i+ E 

n/2+l—drpn—d- 



3 + d-l 
d-l 



j + d-2 
d- 1 



T 



n—d—l 



3 + d-l 
d-l 



q 



j + d-2 
d-l 



(I 



n/2+l—d—jrpn—d—j—l 



By substitution n — d — j — I = i, (14.71) equals 

~n-d-2 



N, 



r.q 



1 _|_ q r-n/2 



E 

i=0 



n — i — 2 
d- 1 
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The formulas (j4"3]h (jOl) and g^D give 

"n-d-2 



iV r ,, 



1 + <f -™/ 2 



E 

i=0 



I n/2+l—drpn—d—l _|_ 



n — i — 3 
d- 1 



n — i — 2 
d-l 

n— 4 

E 

i=d-2 



i+2—n/2rpi _|_ n/2+l—drpn—d—l 



n — i — 3 
d-l 



cZ-3 



z + 2 
d-l 
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i + 1 
d- 1 



j+2-n 



/2yi 



(4.6) 



(4.7) 



(4- 



(4.9) 



We make Fi(T) by gathering positive terms in (14. 9 p and F%(T) from negative ones. | 

Theorem 11.61 claims that all the roots of P r>9 (T) above lie on the circle \T\ = l/y/q if q > 4. 
This is proved in several steps. We consider "normalized" zeta polynomial P r . j(? (T/ v /g). Then 
the Riemann hypothesis is equivalent to the fact that all the roots of P r , q {T '/ \/q) lie on the unit 
circle. On the other hand, P r , q {T '/ \/q) is self-reciprocal, which is the result of the functional 
equation (II. 3p (E^ =0 ajT* is called self- reciprocal if = a v -i for all i). Moreover, if q > 4, 
P r , q (T '/ \/q) turns out to be of the form 



P r , q ( T /Vq) =a + ai T + --- + a k T k + a k T m ~ k + a k ^T m ~ k+1 + ■■■ + a T 



10 



with m > 2k and ao > a% > ■ ■ ■ > a k > 0. We can prove that all the roots of a self-reciprocal 
polynomial of this form lie on the unit circle using several results of classical function theory 
(an analogue of the Enestrom-Kakeya theorem, see Theorem 15. ip . We state the proof in the 
next two sections. 

Remark. We can also prove directly that F\(T/y/q), F 2 (T/yfq) and P T ^ q (T / y/q) are self- 
reciprocal, using the expressions in Theorem 14.51 

5 An analogue of the Enestrom-Kakeya theorem 

In this section, we prove Theorem 1 1.71 This is a self-reciprocal analogue of the following theorem 
and our proof of Theorem 11.71 is based on it: 

Theorem 5.1 (Enestrom-Kakeya) Let f(T) = a + aiT + - ■ - + a k T k satisfy a > a\ > ■ ■ ■ > 
af. > 0. Then f(T) has no roots in \T\ < 1. 

Proof. Marden [TOl p. 151, Exercise 4]. | 

Now, suppose a self-reciprocal polynomial 

f(T) = a + a 1 T + --- + a k T k + a k T m - k + a k . 1 T m - k+1 + --- + a T m (m > 2k) (5.1) 

satisfies a > a x > • • • > a k > 0. We write f(T) as a sum of two polynomials P(T) and Q{T): 

P(T) := a k T m - k + a k ^T m - k+1 + --- + a T m , 

Q(T) := a + a{T + • • • + a k T k , (5.2) 

so f{T) = P(T) + Q{T). Then, by the assumption a > a\ > ■ ■ ■ > a k > 0, we can see from 
Theorem 15.11 that Q{T) has no roots in \T\ < 1. We apply Rouche's theorem to f{T). For 
simplicity, we state it in a restricted form: 

Theorem 5.2 Let C be a circle in C, D be the inside of C . Suppose functions P{T) and Q(T) 
are holomorphic in C U D and \P(T)\ < \Q{T)\ on C . Then Q(T) and P(T) + Q(T) have the 
same number of zeros in D. 

Proof. Ahlfors [Tj p. 153, Corollary]. See also Lehmer [HI Lemma 3]. I 

For our polynomials P{T) and Q{T), we can prove the following: 

Theorem 5.3 We have \P(T)\ < \Q(T)\ on \T\ = r for any r with < r < 1. 

By this theorem, we can see that Q(T) and f(T) = P(T) + Q(T) have the same number of 
roots in \T\ < r. By Theorem 15.11 again. f(T) has no roots in \T\ < r. Since r is arbitrary 
in < r < 1, we can verify that f(T) has no roots in \T\ < 1. Now recall that f(T) is 
self-reciprocal. We have 

T m f (Vj=f(T). 

From this formula, we see that there is a one-to-one correspondence between a root in \T\ < 1 
and that in |T| > 1. We can conclude that f(T) has no roots also in |T| > 1, and all the roots 
of f(T) lie on \T\ = 1. Hence we get Theorem 11.71 

Proof of Theorem 15.31 

First we need the following: 
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Lemma 5.4 (Lagrange's identity) For any Ai,Bi e C, we have 

i=0 i=0 i=0 0<i<j<k 

Proof. Ahlfors [1, p. 9, Exercise 5]. | 

Using this, we can prove the following: 

Lemma 5.5 For P(T) and Q{T) in 115. 2\) . we have 

\P(T)\ = \Q(T)\ 

on \T\ = 1. 

Proof. By letting Ai = and Bi = T % in Lemma [5 A\ we get 

\Q{T)\ 2 = (fc + l)(ajj + ... + a 2)_ ^ l^-a^'-f (5.3) 

0<i<j<fc 

since |T| = 1. As to P(T), noting that |P(T)| = \a k + a-k-iT + ■ • • + doT k \ on |T| = 1, we have 
\P(T)\ 2 = (k + l)(a 2 + --- + a 2 k )- W-i~a k ^T^\ 2 (5.4) 

0<i<j<k 

by letting = a^-i and -Bj = T % in Lemma [5.41 By change of suffices in the sum in (15.41) . we 
have 

\a k - i -a k - j T j - i \ 2 = \ai'-aj'T i '- j '\ 2 = K' " ^"f . (5.5) 

0<i<j<k 0<j'<i'<k 0<i<j<k 

We compare the term for in (15.31) and (|5.5p : 

|aj — OjT^ l \ 2 — \dj — aiT^ % \ 2 = (cij — ajT^~ l ){a,i — ajT J ) — (dj — aiT^~ % ){aj — aiT 3 % ~) 

= 

since |T| = 1. We see that the sums in the right hand sides in (I5.3P and (15 ,4p are the same, 
and we obtain \P(T) \ = \Q(T)\ on |T| = 1. | 

The proof of Theorem 15.31 is completed by invoking the following well-known result: 

Theorem 5.6 (The maximum principle) Let g(T) be holomorphic and nonconstant in a 
bounded (open) region D C C and continuous in D (the closure of D). Then \g(T)\ has its 
maximum M on D — D and we have 

\g(T)\<M 

in D. 

Proof. Ahlfors [1, p.134]. | 

We apply Theorem EH to g(T) := P(T)/Q(T) and_L> := {T E C ; \T\ < 1}. Clearly g(T) is 
meromorphic and nonconstant. It has no pole in D by Theorem 15.11 Moreover, from Lemma 
15.51 |g(T)| = 1 on the boundary of D. Therefore ^(T 1 )! < 1 in D by Theorem 15.61 and we get 
Theorem 15.31 
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6 Proof of Theorem 11.6 



In this section, we prove Theorem 11.61 We have from Theorem 14.5 

p " \S) = TV^n { F > (J) - (^) 



where 



% — % — d 3 

2 — 2 — — 2 

Note that n — d — 1 < d — 3 if r > 3 and g > 4. So there is no term of the same degree from two 
summations in F\(T) of Theorem 14.51 Moreover, P r , q {T '/ ' yfq) is self-reciprocal. It follows from 
the functional equation (II. 3p . but we can verify it directly by showing Fi(Tj \fq) and F 2 (T/y/q) 
are self-reciprocal. Hence we can assume P r ,q(T '/ \fq) is of the form (15.11) . Let 

P r , q {T/y/q) =a + ai T + --- + a^^T^ 1 + a n ^ l T d - 3 + ■■■ + a T n -\ (6.3) 



Lemma 6.1 If r > 3 and q > 4, we have 

a n-d-2 > a n-d-l > 0. 

Proof. Recall d = q r ~ l . Using the expressions (16. ip and (16.21) . we have 

a n - d - 2 = q 2 -^(q r - 2 - 1) 

and 

an-,-! = q^ d)/2 . 

Therefore, a„_d-i > and 

a„_,_ 2 - a n - d - x = g (3 - d)/2 {vW~ 2 " 1) " 1}- 
The last expression is positive if r > 3 and q > 4. | 
Lemma 6.2 If r > 3 and q > 4, we have 

®>i > a i+l 

for < i < n — d — 3. 
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Proof. Using the expressions (16.11) and (16. 2p . we have 



a i+1 = ^-^^-^-^ 



From these formulas, we have 



g (i-n)/2+2 I n * A I ^ _ d _ . _ _ ^ _ . _ 2 ^ _ 



d-l 

= (n — z — 2)(n — i — 3) — (n — i — 3)(n — d — i — l)(q + ^/g) + (n — d — z — l)(n — c? — i — 2)q x /q. 

It suffices to show that the right hand side is positive. It is a quadratic function of the parameter 
i, so we denote it by g(i) = ai 2 + hi + c. We can show that a, 6, c > if g > 4. Indeed, first we 
have 

a = q^q + 1 - (q + y/q) = (y/q ~ l)(q ~ 1) > 
if q > 2. As to b, recall n = (q r — l)/(g — 1) and d = q r ~ 1 . We have 

V^(g - 1)6 = q r (^q - l)(q - 1) + q 3 / 2 (3q 3 / 2 - Aq - 5^ + 7) + ^(2^ - 3) 

(such calculation can be easily done with the help of some expression manipulation program). 
As above, q r (\/q — l)(q — 1) > if q > 2 and 2^/q — 3 > if q > 3. We can easily show 
3 g 3/2 _ Aq _ 5 ^ + 7>0ifg>4 (show that (3g 3 / 2 -Aq- 5^/q + 7) \ q=A > and (3g 3 / 2 - Aq - 
5y/q + 7)' > in q > A). Therefore b > if q > A. 
We can similarly show c > 0. Because 

^(g-l) 2 c = q r{ q ( q ^-2q-2q 1 / 2 + A) + (q 1 / 2 -2)} 
+ g 5 / 2 (2g 3 / 2 -3g-4g 1 / 2 + 7) 
+ g 1/2 (g 2 + 2g 3/2 -7g + 2), 

and we can show that g 3 / 2 — 2g — 2g : / 2 + 4 > if q > 4, and that all other functions in the 
parentheses ( ) are positive in q > A. 

It follows that g(i) > if i > and d{ — a i+ i > 0. | 

We can conclude from Lemmas 16.11 and 16.21 that 

a > ai > ■ ■ ■ > a„_ d _i > (6.4) 

for the coefficients in (16.31) . The assumption of Theorem II .71 is satisfied and the proof of Theorem 
11.61 is completed. 

Remark. (1) If we write W^Ham(r,g) { x -> V) m the form (11.11) . n = (q r — l)/(? — 1) and d = 3. 
Thus we have found infinitely many invariant polynomials satisfying the Riemann hypothesis 
for a small d. 

(2) When q = 2, 3, the coefficients of P r>q (T/ */g) does not satisfy (16. 4p as the following examples 
show. So, in these cases, we cannot prove the Riemann hypothesis in a method described so 
far, but numerical experiments imply that it is very plausible that the Riemann hypothesis is 
true also for q = 2, 3. 
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Example 6.3 (i) Let r = 3, q = 2. Then 



7 7 

WHam(3,2)(x, y) = X 7 + y=xY + 7xY + -= V \ 

l + v 2 l + v 2 



F 1 (T)-2F 2 (T) = -L + (1 + V2)T + (2 + v^T 2 + 2T 3 , 

v2 



v^y " Vv^y V v^y V v^y v 7 ^ 

l ._ (T+1) ( T -=l±i) ( T -- 1 -' 



V2 V y V Vz 

Hence W / Ham(3,2)(^, y) satisfies the Riemann hypothesis. 

(ii) Let r = 4, g = 2. Then P^T/y/q) is of degree 11. We normalize P 4j2 (T/y/q) with a 

suitable constant C as C P^ 2 (T / y/q) = 1 + aiTH . Then we can approximate the coefficients 

as follows: 

a = 1 a 3 m 1.028518954 

ai « 1.414213562 a 4 « 0.606060606 
a 2 « 1.363636363 a 5 ~ 0.317735799 

and a 6 , • • • , an are the same as above in the reverse order. We have a < ai > a 2 > a 3 > 
a 4 > a$ > 0, but according to the numerical experiment, W / Ham(4,2)(^, y) seems to satisfy the 
Riemann hypothesis. 

(iii) Let r = 5, q = 2. Then P 5;2 (T/^q) is of degree 27. Choose C as CP 5;2 (T/y/q) = 
1 + a\T + • ■ ■. Then the coefficients are approximated as 



a 


= 1 


a 7 


« 0.2623468638 


ai 


« 1.414213562 


a 8 


« 0.1391304348 


a 2 


w 1.444444444 


a 9 


« 0.0655867159 


a 3 


« 1.257078722 


a 10 


« 0.0268497330 


a 4 


« 0.977777778 


an 


« 0.0092051531 


a 5 


« 0.691393297 


Ol2 


« 0.0024887453 


a 6 


« 0.446376812 


ai3 


« 0.0005216551 



and a 14, • • • , a 2 7 are the same as above in the reverse order. In this case we have a < a x < 
a 2 > a 3 > • • • > a 13 > 0. The Riemann hypothesis seems to be true. 

Example 6.4 Let r = 3, q = 3. Then P 3t3 {T/^q) is of degree 9. Choose C as CP 3 , 3 (T/ v /g) = 
l + aiT+- ... Then ai « 1.039230485, a 2 = 0.6, a 3 « 0.1732050808, a 4 = a 5 = 0, and a 6 , • --,09 
are the same but in the reverse order. The Riemann hypothesis seems to be true. In many 
other Ham(r, 3) with r > 4, we can observe that the coefficient of T in P r ^{T / yfq) is greater 
than the constant term. 
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7 The Golay codes 

In this section we consider the case where C = Q23, the binary [23, 12, 7] Golay code or C = Qn, 
the ternary [11,6,5] Golay code. We have 

Wg 23 (x, y) = x 23 + 253a; 1 V + 506x 15 y 8 + 1288x 12 y u + 1288x n y 12 + 5O62V 5 

+253x 7 y w + y 23 , (7.1) 
Wg u {x,y) = x n + 132xV + 132xV + 330xV + 110xV + 2% n (7.2) 

(see pH p.94]) or 

Wg 23 ± (x, y) = x 23 + 506x 15 y 8 + I288x u y 12 + 253xV 6 , (7.3) 
W Qll ±(x,y) = x 11 + 132xV + HOxV- (7.4) 

The case C = Q\\ is quite easy: 

Proposition 7.1 The zeta polynomial Pg 11 (T) of the invariant polynomial Wg u (x,y) is given 
by 

/3 — 1 

PQn(T) = ^^(v / 3T+l)(3T 2 + 3T+ 1). 
All the roots of Pg n (T) lie on the circle |T| = 1/V3. 

Proof. The explicit form of Pg n (T) can be obtained by computer calculation. The latter 
statement is obvious. | 

Next we consider C = Q<zz- By computer calculation we get the following: 

Proposition 7.2 Let Pg 23 (T) be the zeta polynomial of the invariant polynomial Wg 23 (x,y). 
Then 

25194 ~ f T 



2-V2 



P&3 (^) = 13 ( 1 + Tll ) + ( 13 + 39 ^2)( T + Tl0 ) 



+ (130 + 39^) (T 2 + T 9 ) + (130 + 156^) (T 3 + T 8 ) 

+ 591 ± 3^ r + r7) + 591 + 459V2 (r5 + n ^ 

The polynomial Pg 23 (T/\/2) does not satisfy the assumption of Theorem 11.71 We would like 
to verify the Riemann hypothesis for Wg 23 (x, y) as theoretically as possible. Our method is 
influenced by that of Duursma [HI Section 5] and [3 Section 5]. 

Let f(x) be a real polynomial of degree n. Then f*(T) = T n f((T + T _1 )/2) is a self- 
reciprocal polynomial of degree 2n. If T = e l9 , then (T + T _1 )/2 = cos#, so the behavior of 
f*(T) on the unit circle can be captured by the behavior of f(x) in the interval [—1, 1]. We 
denote this mapping by p: 

l> :!-f. (7.6) 

We would like to pull back f(x) from a given self-reciprocal polynomial f*(T). It turns out 
that the inverse mapping p" 1 always exists. To clarify this, we introduce two linear spaces of 
polynomials: 

V n := {a + a\X + • • • + a n x n ; a,- G R}, 



3 

W n := {b + b{T + • ■ ■ + b n T n + 6 n _ 1 T n+1 + ■ ■ • b T 2n ; b 3 G R}. 
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The operations in the vector spaces are the same as ordinary summation of polynomials and 
multiplication by real numbers. 

Lemma 7.3 The mapping p : V n — > W n defined by l[7.6\ ) is a linear isomorphism. 

Proof. Clearly p is linear. For the later use, we describe the matrix A p of p with respect to 
the bases 

V n = R[l,X,---,X n ], 

W n = R[T n ,T n ~ 1 + T n+1 ,---l+T 2n ]. 
Because p maps 1, x, x 2 , ■ ■ ■ as 
1 i-> T n , 



2 

2\ (2\ 



x 2 i-> -jr(T n ~ 2 + T n+2 ) + -jrT n , 

x 3 i-> ^(T^ + T^ + Q-iT^ + T^ 1 ), 

x 4 i-> ^.(T n ' i + T n+4 ) + ^-(T n - 2 + T n+2 ) + ^-T n , 



we have for even n, 



1 

2 -Ui 



G) 



2-0 






2- g) 


2" 3 Q 



2- 4 G) • 










2" 4 (t) " 


^ \n/2-l) 








2" 4 ffl " 


n—n( n \ 







where all the elements in the lower triangular part are zero. If n is odd, the (n + l)-th col- 
umn is replaced by '[0, 2" n ( (n _ n 1)/2 ), 0, 2- n (, J, • • • , 2" n (™)]. By this expression, det A p = 
2 -n(n+i)/2 _^ o, so A p is regular. | 

We sketch the proof that all the roots of Pg 23 (T/ y/2) lie on the unit circle. The calculation is 
done with the help of a computer. Let F(T) = (25194/(2- V2))Pg 2S (T 2 / ^2). Then degF = 22 
and F(T) is self-reciprocal. Construct A p for n — 11 and map F(T) by p~ l using Ap' 1 . Then 
we have 



26624X 11 + 



-66560 + 19968^) 



x 



(74880 - ZmmV^x 7 + (-45760 + 29952v / 2)x 5 



+ (14324 - 9984 v / 2> 3 + (-1754 + 1239V2>. 

We can also verify (p- 1 F)(-l) < 0, (p- 1 F)(-0.8) > 0, (p- 1 F)(-0.6) < 0, (p- 1 F)(-0.4) > 0, 
(p" 1 F)(-0.2) < 0, (p" 1 F)(-0.1) > 0. It follows that (p- x F)(x) has at least five roots in the 
interval (—1, 0). It has the same number of roots in (0, 1) since it is an odd function of x, and 
(p _1 .F)(0) = 0. Thus all the roots of (p _1 F)(x) are distinct and lie in (—1,1). Let T = e td . 
Then x = cos 8. While x moves from 1 to — 1, T 2 goes once around the unit circle. Therefore 
Pg 23 (T 2 /i/2) assumes zero exactly 11 times on |T| = 1. 



17 



8 Appendix — an elementary proof of existence of P(T) 

Existence and uniqueness of the zeta polynomial P(T) for a linear code was first established in 
[H Section 9], but a detailed proof is not given. Here we give an alternative, elementary proof, 
including the case W(x, y) G C[x,y}. 

Suppose W(x, y) G C[x, y] is a polynomial of the form (11. II) . First note that 

/CO := 7 77 -T) + xT) n 

J\ ) (l_T)(l-gT) VyV ; ; 

= (1 + T + T 2 + ■ ■ -)(1 + qT + q 2 T 2 + • • -)((x - y)T + y) n 

n 

= (1 + Cl T + c 2 T 2 + • ■ ■) C) {X ~ !lY!l " } 



for some c, G N. Expanding the last formula, we find for some integers b 



y ■ 



the constant term = y n , 

the coefficient of T = nxy n ~ l + (cx — ^)y n , 



the coefficient of T* = biox l y n % + bnx % 1 y n 1+1 + • ■ ■ + &ii£/ n , 

the coefficient of T n ~ d = b n _ d>0 x n - d y d + b n _ d>1 x n - d - l y d+l + ■■■ + b n ^ n _ d y n . 

Let ao, ai, • • •, a n -d G C and we form a function F(T) := (a + aiT + • • • + a n -dT n ~ d ) f (T) . 
Then the coefficient of T n ~ d of F(T) is 

+a n ^ d _ 1 {nxy n ~ 1 + (ci - n)y n } 



+a,{& i0 xV- J + bax'-V^ 1 + ■■■ + buy n } 

+a {b n „ dfi x n ~ d y d + b n ^ d>1 x n - d -y +1 + ■■■ + b n ^ d y n }. (8.1) 

On the other hand, since (W(x, y) — x n )/(q — 1) = (A d x n ~ d y d + ■ • • + A n y n )/(q — 1), we can 
determine ao, ai, • • •, a n -d so that (18.11) coincides with (W(x,y) — x n )/(q — 1) (the system 
of linear equations for determining ao, a 1; • • -, a n _ d has a regular coefficient matrix). So we 
can always determine the zeta polynomial P(T) from a given W(x, y) uniquely as P{T) = 
a + a{T + • ■ ■ + a n _ d T n ~ d . | 
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